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We present a finite temperature calculation of the transverse spin-diffusion coef- 
ficient, D±, in a dilute degenerate Fermi gas in the presence of a small external mag- 
netic field, H. While the longitudinal diffusion coefficient displays the conventional 
low-temperature Fermi- liquid behavior, oc T~^, the corresponding results for D± 
show three separate regimes: (a) D± ~ H^^ for T <^ H; (b) D± ^ T^'^, D_\_/D\\ ^ 1 
for T ^ H and large spin-rotation parameter ^ ^ 1, and (c) D± = oc 
for T ^ H and ^ ^ 1. Our results are qualitatively consistent with the available 
experimental data in weakly spin-polarized ^He and ^He — ^ He mixtures. 
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The unusual features of spin dynamics in spin-polarized quantum systems have been 
intensively studied since the pioneering paper of Leggett and Rice |1| on spin diffusion in 
normal liquid ^He. The main effect arises from the observation that the presence of a molec- 
ular field (induced by the applied magnetic field) leads to the an additional precession of the 
spin-current which in steady state acquires a component perpendicular to the magnetization 
gradient; through the continuity equation this results in an anomalous reactive component 
(damped spin-wave) to spin transport. This "spin-rotation" effect is also present in the case 
of spin-polarized Boltzmann gases |^. 

From a microscopic point of view, it was expected that a number of qualitatively novel 
phenomena might arise in the case of sufficiently high polarizations 0. A natural suggestion, 
made by Meyerovich 0], was that low-temperature spin- diffusion becomes highly anisotropic 
for finite polarizations. More precisely, processes involving spin-flips make use of the phase- 
space volume between the two distinct Fermi surfaces (for "up" and "down" spins) leading 
to a finite scattering rate in the limit of T ^ 0. This is in contrast with processes involving 
scattering in the vicinity of each of the Fermi surfaces which are subject to the phase-space 
restrictions of unpolarized Fermi liquids and are thus characterized by the conventional 
Fermi-liquid behavior of scattering rates, oc 0. Meyerovich's suggestion was recently 
supported by measurements of the transverse spin-diffusion coefficient in weakly polarized 
liquid '^He 0. Theoretically, was calculated in the dilute gas limit at T = [^|[ and only 
a tentative estimate based on a variational solution of the Boltzmann equation is available 
for T 0. 

In this letter we present the analytical calculation of the finite temperature behavior of 
D± for a dilute, weakly polarized Fermi gas, in the s-wave approximation by solving the ap- 
propriate kinetic equations exactly. To the extent to which in the limit of small polarizations 
strong interactions only lead to constant renormalizations of the weakly interacting result 
we expect that our findings should also apply to the case of weakly polarized ^He. Indeed, 
as explained below, the detailed low {H, T) behavior of our expression for is consistent 
with the small systematic deviations (lying within the error-bars) in the results of Ref. |^ 
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from previous theoretical estimates. We note that the exact solution of the kinetic equations 
does not merely lead to quantitative renormalizations of the transport coefficients; rather, 
it brings out new qualitative effects. In particular, due to the existence of two dimensionless 
parameters, the low-temperature behavior of transverse diffusion displays two crossovers. 



the first at T ~ and the other at T oc ^/Hep {ep is the Fermi energy). Our finding, not 
appreciated in previous work on the subject, is that spin diffusion becomes isotropic (i.e., 
D± = D\\) only for T above the second crossover. 

Our starting point is the kinetic equation for the Wigner transform of the transverse 
component of the density matrix (in the frame rotating at the bare Larmor frequency), 
nii{t, f,p), where t, r are the center of mass time and space coordinates and p is the relative 
momentum vector: 
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Here nt'^ = [exp((ep =F H/2 — tF)/T) + are the equilibrium distribution functions for 
up- and down-spin particles, respectively, a is the s-wave scattering length, and ep = p'^/2m. 
Throughout this paper we use energy units for the field H and we set h = 1. The two 
terms in the r.h.s. of equation (|l|) correspond to direct and exchange two-particle scattering 
processes, described by the finite-temperature many-body T-matrix, 
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where {g^, g) is the 4- momentum of the center of mass, the incoming particles have momenta 
g±Pi, and the outgoing - g±p2. 

In the kinetic equation (|I]), the functions A{p,p') and B{p,p') contain all effects of 



second order in app, Pf = y/2meF- In the case of low temperature and small polarization, 
T,H <^ eF, these take the form. 
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A more detailed derivation of the kinetic equations will be given in a longer publication. 

To compute the transverse spin-diffusion coefficient we will solve the kinetic equations 
(HI) and extract the steady state transverse spin-current driven by a constant magnetiza- 
tion gradient. For simplicity we will consider the case in which Mj_{r) = Mj_{x) with 
dMj_{x)/dx = const. We will also restrict ourselves to the low temperature, weakly polar- 
ized limit, T, H <^ ep, in which case a complete analytical solution is possible in the s-wave 
approximation. 

With these assumptions the solution of the kinetic equation (|lD can be taken to be of the 
form, n^i{t,f,p) = g{t,x,p) + f{t,x,p) cosip, where p = \p\ and cosip = p- x/p. Performing 
explicit integrations in equation (|^) then leads to the two coupled kinetic equations, 
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In (H-^ Irel and Isr represent the relaxational and spin-rotation parts of the "collision 
integral" , 



Irel[Wip)] = W{p) A{p,p 



(27r)3 



A{p',p)W{p 



\27r)3 



(7) 



Isr[W{p}] = -W{p) J Bip,p'){ni,-ni,)0^ + {ni-ni) J B{p,p')W{p 
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As can be shown by analysing the eigenvalue spectrum of the relaxational part of the collision 
integral Ireii[f] the spin current decays to its steady state value in a microscopic time scale, 
T±, beyond which the time derivative terms in equations can be omitted. In addition, 

since dM±{x)/dx = const or equivalently dg/dx is independent of x, equation implies 
df /dx = 0. The solution of can then be shown to be 

go{x,p) = G{x){nl-ni) . (9) 

in which case the equation for f{p) can be written as 

VF^{nl-4) - = -Irelllfip)] ' ^Isrl[f{p)]. (lO) 

Here, vp is the Fermi velocity and fi^^-* = (87ra/m)M|| ^ 2apFH/'K represents the leading cor- 
rection to the precession frequency of the spin current, = J± — iJj_ = /o°° f{p)p^dp/QTT^m 
(Af|| is the longitudinal magnetization induced by the field H). Note that no such correc- 
tion appears in the equation for the transverse magnetization itself, Mj = Mf^ — iM\_ = 

g{p)p'^dp/ An'^ , i.e., there is no renormalization of Larmor precession, as expected from 
rotational invariance of the inter-particle interactions. 

To identify the diffusion coefficient we want to match the conventional form of the 

macroscopic constitutive relation, 

to the one implied by the microscopic equation (p^O]). Here, is the transverse spin diffusion 
coefficient while ^ is referred to as the "spin-rotation" parameter. The non-zero imaginary 
part of the denominator in (|Tl|) reflects the fact that, due to the spin- rotation effect, the spin 
current is not parallel (in spin-space) to the driving magnetization gradient. We recall that 
in the relaxation-time approximation |1|] ~ ^^^''Treiax, where Treiax is the corresponding 
relaxation time. We note from the outset that, except in the high field regime, H 3> 
T, the relaxation-time approximation breaks down and the above parametrization of ^ is 
inapplicable. We are now in position to compute D± and ^ as functions of T and H. As we 
show below, the calculation is tractable analytically in both "high"- and "low" -field limits: 



High-Field Behavior, H/T S> {apF)iT /^f)'- In this case the relaxational term in (|T0]) is 
small, Ireii[f] ^ ^^^^Z, and the current is almost perpendicular to the magnetization 
gradient. In this limit the solution of equation (|10|), fnip), can be obtained iteratively: 

Mp) = -^^^(4-4) - • (12) 

Upon integration this leads to 



(13) 



where J±[Vr(p)] = / Vr(p)p3rfj9/127r2m. In (|g) and (0) we identify J±[/srl[4-^^i-]]/J± hi- 
nt] = — with the second order correction to the spin current precession frequency, and 
J A^reixW^p — "n^^l J ^\n',^ — nt] = td{H,T) with the "diffusion time". We note that the 
"diffusion time", td, is not the relaxation time of the spin current usually involved in the 
phenomenological discussions: the latter describes the relaxation of the distribution function 
towards the steady state solution, foip) and would arise as the eigenvalue of the relaxational 
part of the collision operator (see below). With these definitions, (|TT|) and (13) immediately 
yield. 



The first thing to stress is that the behaviour of D± in ([1^) differs significantly from 
that of the longitudinal spin diffusion coefficient [jlT| , D\\ ^ (f|,/87rma^T^)C(— 1/3), where 



C(— 1/3) ~ 0.843 is the Brooker-Sykes coefficient [0. The origin of this effect is the different 
phase space restrictions associated with the scattering processes leading to transverse and 
longitudinal spin diffusion in the high- field limit. While collisions responsible for longitudinal 
spin diffusion are restricted to energies within fc^T of each of the ("up" or "down") Fermi 
surfaces, those leading to transverse spin diffusion involve spin-flips which can also take 
advantage of the full region enclosed between the two surfaces. Thus, the expression for D± 
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in (|TJ) involves the phase-space for scattering generated both by the magnetic field as well 
as by the thermal smearing of the individual Fermi surfaces. As a result, the transverse 
diffusion coefficient in (|1^) remains finite at T = as was predicted by Meyerovich |^ and 
recently confirmed experimentally 0. 

It should be noted that, even ioi T ^ H (but still in the high field limit, H ^ 
{apF)T'^ /ep) the transverse spin diffusion coefficient D± ^ {3vp/32TTma'^T^) ^ 0.890D|j 
still differs from the longitudinal one. This is contrary to the commonly implied belief |]7|JT5 



that the crossover to the isotropic regime of spin diffusion occurs for T ~ iJ. As we will see 
below D± and D\\ become equal only at much lower fields, H ~ {app)T'^ /ep- 

Crossover and Low-Field Behavior, H/T ~ {apF)T /ep'- In this region, Ireii[f{p)\ ~ 
qS^) f[j)^ ^ Isri[f{p)] and, therefore, the spin- rotation terms in the collision integral may 
be omitted. Furthermore, one may also set = in evaluating the functional Ireii- In 
terms of reduced variables, rj{j)) = (cp — ep)/T (and f{rj{j))) = f{p)), h = H/2T ^ 1, and 
Fq = 2a^m/37i, the steady state equation (0) then becomes, 

-2ro/_> + a)/(.)(-^ + -^)... (16) 



Our analysis makes use of the methods of Reference |]T2|] and begins by transforming (|TBp 
into the differential equation, 

n 9„,,, 2 , n , nhvp dG coskh . 
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for the function F{k) = /^^ e^^V(^?) cosh(r//2)dr/. Here 7^ = 1 - in'-l)/37r'^ToT'^ ^ 1 - 
2iH€p/TT^T^app. 

The solutions to ([17|) may be expressed in terms of Gegenbauer polynomials ||l3l with a 
complex index, 

F{k) = J2 FnMk) ,Mk) = (cosh7rfc)-^Cn^^(tanh7rfc). (18) 

ra=0 
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Note, that since the functions f{ri) and F{k) are even and /c) = (— l)"'0n(/c), only even 
n terms contribute to the sum in ([18D . 

The transverse spin current can then be evaluated as 
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This allows us to write the final expression, 



(22) 



47rma2T2;^o| + (7 + n)(7 + n + l)' ^ ' 

from which we can extract and the spin-rotation parameter ^. The low field limit, 
H « {apF)T'^/€F corresponds to 7 —> 1 in which case ( pSf ) yields, D± = D\\ = 
{vp/8'n"ma'^T'^)C{—l/3). In this limit the spin- rotation parameter is given by ^ ~ 
0.139HeF/T'^{app) << 1, refiecting the fact that the spin current becomes parallel (in spin 
space) to the driving magnetization gradient. The detailed behavior of D± as a function of 
T/H, which clearly displays the two crossovers discussed above, is shown in Figure 1. 

Although the calculations presented above are restricted to low densities, we expect that 
our qualitative arguments, based on the existence of two independent parameters, ^ and 
H/T, should also hold in the Fermi liquid regime, T,H << tp- The s-wave approximation, 
which made our analytical calculations possible, raises more serious issues, especially con- 
cerning the detailed behavior of D_l in the crossover regions. In particular, we would expect 
that the region in which oc T"^ with D± 7^ D\\ should shrink. 
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Some comments are in order concerning the possible relevance of our findings to experi- 
ment. First we note that the available experimental data in weakly polarized ^He ^ deviate 
systematically from the naive theoretical fit which uses a single adjustable parameter (T^ in 
0) to cover the entire temperature range including both, ^ > 1 and < 1. Much better 
agreement is obtained by restricting the fit to the ^ < 1 region with an overall prefactor 
smaller than the one implied by fitting to the value of D\\ in the low field, high temperature 
regime (we use for the inverse overall prefactor the value A ^ 7.7 x lO^sec/cm^K^ instead 
of 5.8 X lO^sec/cm^K^ as in Reference |^ and for the "anisotropy temperature" ^ 12.5 
mK instead of 16.4 niK). This is consistent with our picture, with D± < D\\ for > 1 
and T » H. In addition, although it appears that the region between the two crossovers 
cannot be clearly identified - most likely due to large Fermi liquid renormalization effects - 
the isotropic limit is indeed reached in the regime ^ < 1 0. In principle, our calculations 
should be more relevant to the measurements in dilute ^He -^He mixtures. Although in 
the available data (for .18% ^He) the crossover to the isotropic limit occurs for ^ ~ 1 with 
H << T, the temperature is not sufficiently far below ep and, moreover, the polarization 
is somewhat high, ~ 25%. Nevertheless, for ^ > 1 D± oc D\\ with the ratio D^/D\\ slightly 
less than unity |T^,ITB|. Also, the measured T dependence of the "spin- rotation" parameter. 



^, near the crossover to the isotropic ( "low- field" ) limit is qualitatively consistent with our 
results in both the data of reference fl^ and those obtained in the degenerate regime of 
more concentrated solutions (2.6% '^He) [jl^ with lower polarization (~ 2%). In both situ- 
ations, the crossover to the isotropic regime can be clearly distinguished. However, in the 
former case the polarization was rather high (~ 25%) and the "low-field" crossover itself 
occurs beyond the degenerate limit. Also, there is a large discrepancy in the magnitude of 
the shift of ^T^ in the latter case (see Figure 2) which can be attributed to Fermi liquid 
renormalizations anticipated in high concentration solutions. To sharpen the identification 
of two crossovers the data of reference |]1^ should be extended to lower temperatures (to 
study the H/T ~ 1 behavior). Quantitative comparisons could be made only in the more 
dilute case of Reference fl^ where lower field and lower temperature experiments should be 



performed. 

We close with a brief comment about the time scale for the relaxation to the steady 
state solution. This relaxation proceeds exponentially with a transverse relaxation time t± 
defined (in the appropriate reference frame) by 

--fip) + id^^^fip) = -IrelAf] - lIsAf] , n^''' ~ n'^'^ . (24) 

We have only solved this equation in the limit H ^ T'^/^f in which case the spin- 
rotation term Igri results in the rapid oscillations of the distribution function f{p) with a 
characteristic frequency ~ fi*^^^ >> 1/t±. Upon averaging these oscillations over time scales 
short compared to t±, ( ^4]) reduces to an eigenvalue problem for the collision operator, Ireii, 
f{p)/Ti_ = Ireii[f{p)]- The result of this calculation, plotted in Figure 3, illustrates two 
important points: first of all, t± is different from the "diffusion time", td. Moreover, its 
T, H dependence cannot be reduced to the "conventional" form, rj^ oc AT'^ + BH^. 
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FIGURES 

FIG. 1. Spin diffusion coefficient in tlie liigh-field region H ^ T. The transverse spin diffusion 
coefficient D± (equation (|l^)), the longitudinal coefficient Dy, and the naive fit + 
for D± {D^^^ = STTVp/Sma'^H'^ being the limiting value of D± at T ^ 0) are represented by the 
solid, dashed and dotted lines, respectively. The inset shows the ratio D±/D^^ near the second 
(low-field) crossover, T = \/Hep/app. 

FIG. 2. The low-field behaviour of the spin-rotation parameter ^ in the concentrated ^He-^He 
mixture. The points represent the experimental data of Reference jl^, the dashed line - the 
fit for the behaviour of ^ below the crossover, and the solid line is our theoretical result for the 
s-wave approximation (with the values of T and ^ scaled by appropriate factors). Note, that the 
latter is not a straight line. On the inset, the theoretical curve for ^ at T ~ Tc = \J Hep j app is 
plotted in a different scale. 

FIG. 3. The relaxation time of the transverse spin current plotted as a function of T/H. The 
solid line corresponds to the exact eigenvalue of collision integral and the dotted line - to the fit 
rj^ = AT'^ + BH^ , which differs from the exact value in the crossover region H ^ T. The dashed 
line represents the diffusion relaxation time td (see equation (0)). 
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